Microscopic approach to entropy production 
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It is a great challenge of nonequilibrium statistical mechanics to calculate entropy production 
within a microscopic theory. In the framework of linear irreversible thermodynamics, we combine the 
Mori-Zwanzig-Forster projection operator technique with the first and second law of thermodynamics 
to obtain microscopic expressions for the entropy production as well as for the transport equations 
of the entropy density and its time correlation function. We further present a microscopic derivation 
of a dissipation functional from which the dissipative dynamics of an extended dynamical density 
functional theory can be obtained in a formally elegant way. 
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I. INTRODUCTION 

One of the central challenges of statistical mechanics 
is to calculate and predict the entropy of a given system 
under well-specified conditions from a microscopic point 
of view. In equilibrium, this problem dates back to Boltz- 
mann and it is by now standard textbook knowledge that 
the entropy can be obtained by a suitable phase-space 
averaging in different ensembles [l[. The same problem, 
however, is much more complicated and in general un- 
solved in a nonequilibrium situation: while the second 
law of thermodynamics predicts a global entropy pro- 
duction in a closed system, it is very difficult - even in 
principle - to calculate the entropy production within 
a microscopic approach that starts from the individual 
interactions of the particles. Linear irreversible thermo- 
dynamics 1 3] provides a framework, where this prob- 
lem can be addressed systematically in a simpler way. 
Moreover, the Mori-Zwanzig-Forster projection operator 
technique (MZFT) [3-flCjj , when supplemented with a se- 
lection of relevant slow variables, can be used to derive 
microscopic expressions for the dynamical evolution of 
these slow variables and their time correlation functions. 
When the MZFT is applied to the entropy density, there 
is, however, a principle obstacle since the entropy density 
does not possess any known microscopic operator. This 
fact is opposed to, for example, the particle density and 
the energy density, which both possess corresponding mi- 
croscopic operators. 

In fact, we have recently applied the MZFT in the 
framework of linear irreversible thermodynamics to de- 
rive a generalized (extended) dynamical density func- 
tional theory (EDDFT), which provides a microscopic 
basis for the dynamical evolution equations of the rele- 
vant variables [llj. The approach was only formulated 
for slow conserved variables such as particle and energy 
density. In this paper, we generalize the EDDFT ap- 
proach of Ref. [Tl| to include the entropy density a(f,t) 
as an additional (non-conserved) variable. Assuming the 
applicability of linear irreversible thermodynamics JlJ-Q , 
we derive transport equations for the entropy density as 



well as for its time correlation function C aa (r,r',t,t'). 
The key idea here is to combine the first and second 
law of thermodynamics with the microscopic expressions 
obtained within the MZFT for the remaining slow vari- 
ables in order to circumvent the obstacle that the entropy 
density does not possess a known microscopic operator. 
Thereby, we provide - at least in principle - a link of en- 
tropy production to microscopic expressions. In doing so, 
we also present a microscopic expression for a dissipation 
functional from which EDDFT [11| and in particular the 
traditional dynamical density functional theory (DDFT) 
[l^ - [Tr| can be derived in an elegant way. An important 
application of this reformulation of EDDFT in terms of 
a dissipation functional is the derivation of phase field 
crystal (PFC) models from EDDFT, in particular those 
which involve orientational degrees of freedom [l7l - [20| . 

Our derivation of the transport equations for the en- 
tropy density and its time correlation function is valid for 
a general set of relevant variables including also variables 
that are not considered in the context of current EDDFT. 
Furthermore, the derived expressions are not restricted 
to the hydrodynamic limit (vanishing wave vector k — > 
and frequency uj — » 0) and thus generalize the corre- 
sponding hydrodynamic equations to larger wave vectors 
and frequencies. 

The paper is organized as follows: after preliminary 
remarks about the considered systems and notation in 
Sec. mi we present a derivation of transport equations 
for the entropy density and its time correlation function 
as well as a microscopic expression for the corresponding 
dissipation functional in Sec. IIIII These expressions for 
general relevant variables are afterwards specialized to 
the variables of EDDFT in Sec. IIV1 Finally, we conclude 
in Sec. El 



II. RELEVANT VARIABLES AND THEIR 
TRANSPORT EQUATIONS 

As usual in the context of EDDFT [ll| , we consider a 
grand-canonical (total) ensemble of systems of N parti- 
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cles with the not explicitly time-dependent Hamiltonian 
H(ft]^ and the set t t of phase-space coordinates. These 
systems can be described by the total probability density 
p(t), which is given by the solution of the Liouville-von 
Neumann equation 

'p=-tp=- X -[H,p] , p(t)=e- lt p(Q) (1) 

with the Liouvillian £, the imaginary unit i, the re- 
duced Planck constant h = /i/(27r), and the commutator 
[X, Y]=XY- YX of X and Y. 

A. Relevant variables 

Following the idea of the MZFT, we select a set T t of 
n <C N independent relevant variables (operators) ai(f, t) 
with i G {1, ...,n}, which is sufficient to describe the 
considered system on a much simpler basis, where the 
huge number of irrelevant microscopic variables is pro- 
jected out. Transport equations for the relevant variables 
are given by the Liouville-von Neumann equations 

hi = jC&i = - [H, &i\ , a t (f, t) = e ct ai (r, 0) . (2) 
h 

The set of relevant variables T t is associated with a rel- 
evant probability density p(t). Together with the grand- 
canonical trace Tr, the relevant probability density can 
be used to define the averaged relevant variables 

ai(r,t)=Tr(p(t)ai(r,0)) . (3) 

We assume a given generalized Hclmholtz free-energy 
functional 

F[a x ,...,a n ] = fd 3 rf(r,t) (4) 

Jig? 

with the generalized Helmholtz free-energy density /(r, t) 
that describes the state of the considered system in terms 
of the averaged relevant variables ai(r,t). The general- 
ized Helmholtz free-energy functional (|4]) can be used to 
define the thermodynamic conjugate^ 

a ^ r ^>- 5ai (f,t) (5) 

of the averaged relevant variables a,i (r , t) . In terms of the 
relevant variables and their thermodynamic conjugates, 



1 For consistency, wc denote all quantities that arc directly asso- 
ciated with the set of phase-space variables ft by a hat * , while 
the hat is omitted for all other quantities. 

2 Static stability of the considered system requires that T is 
bounded from below. 

3 The internal energy density is an exception. Its thermodynamic 
conjugate has to be considered differently (see Sees. Ill Bl and 

irrm . 



the relevant probability density p(t) is here specified as 
the generalized grand-canonical probability density 

with the grand-canonical partition sum the (con- 

stant) inverse thermal energy (3 = l/{k-QT), where /cb 
denotes the Boltzmann constant and T is the absolute 
temperature, and the explicitly time-dependent effective 
Hamiltonian 

n 

H cfi {t) = H-J2 d 3 ra\{f,t) ai (r). (7) 
i=i 

The n — n c + n n relevant variables ai(f,t) have to be 
distinguished into n c conserved variables al (r, t) and n n 
non-conserved variables a"(r, t). The same holds for their 
averages a,(f, t), which are distinguished into conserved 
averaged variables af(r, t) and non-conserved averaged 
variables a"(r, t). 

B. Transport equations 

The dynamics of the n c conserved relevant variables 
can be described by the conservation equations 

k c i(r,t) + VfJ {l) (r,t) =0 (8) 

with the local currents jW (r, t) and i G {1, . . . , n c }. Cor- 
respondingly, the dynamics of the averaged conserved rel- 
evant variables is given by 

^(r,t) + V r -.jt J )(r,t) = 0. (9) 

On the other hand, it is assumed that the time-evolution 
of the n n non-conserved relevant variables can be de- 
scribed by the balance equations 

a?(T*t) + = (io) 

with the local quasi-currents $W (f^ t) and i G 
{l,...,n n }, while the averaged non-conserved relevant 
variables are described by 

<(f,t) + $W(r,t) = 0. (11) 

Both the currents (r, t) and the quasi-currents 
can be decomposed into two different contri- 
butions Hl|: 

J« (r, t) = Jf (r, t) + 4 } (r, t) , (12) 

$®(?,t)=$$(?,t)+$®(?,t). (13) 

These are the reversible contributions denoted by the 
subscript "R" and the dissipative contributions denoted 
by the subscript "D" . While the reversible contributions 
are isentropic, i. e., they are not associated with entropy 
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production, the dissipative (reversible) contributions in- 
crease the total entropy of the system. 

Examples for the general relevant variables a£ (r, t) and 
af(f, t) are a concentration field and a local polarization, 
respectively. There are only two possible relevant vari- 
ables that have to be treated explicitly and cannot be 
taken into account by the general set o, (r, t) of relevant 
variables. These special variables are the conserved gen- 
eralized internal energy density e(r, t) or its Legcndrc 
transforms and the non-conserved entropy density a(r, t) 
(see Sec. lIII A|) F1 Since the internal energy density is con- 
served, the transport equations for i(r,t) and s(r,t) can 
be written as 



i(f,t) + V?-J £ {r,t) = 
e{r 1 t) + ^ F -J £ {r,t) = 



(14) 
(15) 



with the internal energy currents J £ (r , t) and J e (r , t) = 

Jr(t^) + Jh^T't)- The entropy density, in contrast, is 
non-conserved and the transport equations for a(f, t) and 
a(r,t) are 



a(r,t) + $' T (f,i) = 
a(f,t) + <$> a (r,t) = 



(16) 
(17) 



with the entropy quasi-currents (r, t) and (f , t) 



C. Dissipation functionals 

If the considered system is in local thermodynamic 
equilibrium so that the local formulation of the first 
and second law of thermodynamics hold, linear irre- 
versible thermodynamics [lHl] can be applied [2l|. This 
useful framework is especially applicable for all passive 
systems and states the existence of a dissipation func- 
tional £R from which the dissipative currents (r, t) and 

quasi-currents (r, t) of the averaged relevant variables 
ai(f,t) and the entropy density cr(r,t) can be derived. 
The dissipation functional 



d 3 r r(f, t) 



(18) 



describes the total amount of energy that is dissipated 
per time in the considered system. Its integrand, the dis- 
sipation function t(r , t) , is therefore positive semi-definite 



4 In fact, the internal energy density i(r,t) has to be treated sep- 
arately, since e(r, t) is a Legendre transform of the generalized 
Helmholtz free-energy density f(r, t) that can be gradient ex- 
panded in terms of the other relevant variables. The entropy 
density a(r,t), on the other hand, has to be treated separately, 
since it is not independent of the internal energy density e(r, t). 
Both variables are related by the local formulation of the first 
law of thermodynamics. 



as a consequence of the second law of thermodynamics 
It is nonlinear in the thermodynamic 
variables and quadratic in the thermodynamic forces that 
arc defined as 



a/ 



c: 



(19) 



and correspond to the averaged relevant variables (r , t) . 
Again, the internal energy density e(r, t) and the entropy 
density o~(r,t) have to be treated separately. While a 
thermodynamic force corresponding to the internal en- 
ergy density does not exist, the thermodynamic force 
cr" (f , t) corresponding to the entropy density is defined 
further below by Eq. For a given dissipation func- 

tional, the dissipative currents and quasi-currents of the 
averaged relevant variables follow directly by functional 
differentiation with respect to the thermodynamic forces 



r« _ 



~5af 



(20) 



The dissipative energy current Jfj (r , t) and the dissipa- 
tive entropy quasi-current $g (r, t) have to be derived 
differently (sec Sec. IIII Al further below). 

For systems far from thermodynamic equilibrium, lin- 
ear irreversible thermodynamics cannot be applied and a 
dissipation functional is in general not known [l(. This 
is the case for active systems like, for example, lasers, 
amplifiers, and biological systems. Indeed there is a gen- 
eralization of the dissipation functional - the generalized 
Ljapunov functional - to systems that can also be far 
from thermodynamic equilibrium, but a Ljapunov func- 
tional does not exist in general. The criteria for the ex- 
istence of a Ljapunov functional for a particular system 
are characterized by potential conditions [2214251 ] . 



III. 



THE ENTROPY DENSITY AS A 
RELEVANT VARIABLE 



The entropy density o~(r, t) is a special (non-conscrved) 
relevant variable that has not yet been considered in the 
context of EDDFT [lj. Although EDDFT could in prin- 
ciple be generalized to include also non-conserved vari- 
ables, the entropy density cannot be incorporated into 
EDDFT directly. It has instead to be considered sep- 
arately, since there is no corresponding operator <r(f, i) 
known for the entropy density that would be required 
to incorporate the entropy density into EDDFT using 
the MZFT [ll|. To avoid this problem, we assume the 
applicability of the local formulation of the first law of 



5 The condition r(r, t) ^ is a necessary prerequisite for dynami- 
cal stability of the considered system [2 ill . For entirely reversible 
processes, r(r, t) is zero, while it is positive for all other (dissipa- 
tive) processes. 
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thermodynamics for operators and express the unknown 
operator for the entropy density in terms of the known 
operators of the other relevant variables. We thus con- 
struct an expression for the entropy density operator that 
is applicable in the framework of linear irreversible ther- 
modynamics. On this basis, we present the derivation 
of dynamical equations for the entropy density o~(r, t) as 
well as for its time correlation function C aa {r,r',t,t'). 



A. Non-equilibrium dynamics 

Since the entropy density o~(r , t) is not conserved, its 
dynamics has to be described by the balance equation 
(|17p . where the entropy quasi-current 



^(f,t) = V r J a (?,t)-Q°(r,t) 



(21) 



can be decomposed into a divergence term with the en- 
tropy current J a (r, t) = J^(r, t) + Jg(f, t) describing the 
transport of entropy and a source term Q CT (r, t) denoting 
the production of entropyH The entropy production 



CT(r,<) 



2r(r,t) 
T(f,t) 



(22) 



in turn can be expressed in terms of the dissipation func- 
tion r(r, t) and the absolute local temperature T(f,t). 
In order to derive explicit microscopic expressions for 
J a (r,t) and Q a (r 1 t) 1 we make use of the local formu- 
lation of the first law of thermodynamic^ [HI 



de = Tda 



da^ 



da? 



(23) 



with the generalized internal energy density e(r, t), which 
is related to the generalized Helmholtz free-energy den- 
sity f(r,t) by 



df = de - Tda - adT 



(24) 



and additional conserved and non-conserved classical 
thermodynamic variables a^(f,t) and af(f,t), respec- 
tivelylfl Equation (f2"3"| can be rearranged into 



1 • a i -c 



nh 

V 

i 

T 



(25) 



6 The dissipative entropy quasi-current is therefore <&q (r, t) = Vf~ 
J£(r,t)-Q°(r,t). 

7 Einstein's sum convention is assumed in the following. Equation 
11231 makes the internal energy density e(r, t) and the entropy 
density <r(r, t) dependent on each other and is the reason, why 
e(f,t) and <r(r, t) have to be treated separately from the other 
relevant variables ai(r, t) (see Sec. Ill Bl . 

8 Notice that the thermodynamic conjugate of the entropy density 
is the local absolute temperature: a^(f, t) = de(f,t)/dcr(f,t) = 
T(r,t). 



providing a dynamical equation for the entropy density 
in terms of the time derivatives^] of the internal energy 
density and the other relevant variables. Using the trans- 
port equations © and (fTTj) for the general relevant vari- 
ables a1(r,t) and af(f,t), respectively, (fTS"]) for the inter- 
nal energy density e{r 1 t), and (|17p for the entropy density 
a(f, t), Eq. (|25[) can be transformed into the balance equa- 
tion for the entropy density (|17p with the decomposition 
(1211) . the dissipative entropy current 



1 



cb, 

J a (r t\ — — T e - ^i- 
J D\ r i z ) — ji J D rp J D 1 



(26) 



and the entropy production 

Q a (f,t) = —JZ-VfT-^.V^f + ^a? . (27) 

Notice that only the dissipative currents and quasi-cur- 
rents contribute to the entropy production so that the 
following equation must hold: 



(28) 



'R "T «i v f J R T *R 

Combining Eqs. (f2"2"j) and (f2~T)) now leads to the desired 
microscopic expression for the dissipation function: 



(29) 



Since the dissipation function t(r, i) has to be positive 
semi-definite, Eq. (|29p implies that the dissipative cur- 
rents and quasi-currents are functions of the thermody- 
namic forces (fT9"|). In fact, these functions are assumed 
to be linear and homogeneous [2l[. With this assump- 
tion, is is obvious that the dissipative currents Jq (r , t) 
and quasi-currents $^(r, t) corresponding to the con- 
served relevant variables aj (r , t) and the non-conserved 
relevant variables a?(r, t), respectively, can be derived 
from a given dissipation functional *H using Eqs. (|20[) . 
Analogously, the dissipative entropy current Jg (r, t) can 
be derived: 



o = -V r -T , 



J 5 - W ■ < 30 > 



The dissipative energy current Jg(r, <) follows then from 
Eq. HMD- 
Equation H7]) together with Eqs. j2T| and ([26 ]) - ([29" ]) 
embody a microscopic expression for the dynamical equa- 
tion of the entropy density. This microscopic equation 
expresses the entropy quasi-current and the dissipation 
functional in terms of the (microscopic) currents and 
quasi-currents of the other relevant variables and consti- 
tutes the first main result of this paper. It is noteworthy 
that the equations in this subsection closely resemble the 
corresponding hydrodynamic equations |21j . 



9 Throughout the whole paper, a dot ' denotes a derivation with 
respect to time t. This applies also to time derivatives of corre- 
lation functions with two time variables t and t' further below, 
where the dot always means a derivation with respect to t and 
not to t' . 
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B. Equilibrium correlations 



tion function: 



We now consider the dynamics of the entropy time 
correlation function C (7(7 (f,f',t,t'). In the linear regime 
near equilibrium, the equilibrium time correlation func- 
tion (Kubo function [8|) of two variables X(r,t) and 
Y{r,t) is defined as @,[n| 



C XY (r,r',t,t') = (AX^(r,t)\AY c( i(r',t')) c 



(31) 



Here, the letters "eq" denote equilibrium quantities, 
(•|-) cq is Mori's scalar product and AX eq (f, t) = 

X(r,t) - X i(f) with A cc >(r) = {X(r,t)) cq are the equi- 
librium fluctuations of the variable X(r,t). Since they 
are associated with equilibrium fluctuations, the time 
correlation functions (f3~Tj) are translationally invariant 
with respect to time: Cxy (r, r', t, t 1 ) = Cxy (r, r', t — t'). 
For the derivation of the dynamics of C a(T (f,r',t,t'), we 
need an expression for the operator a(f, t) of the entropy 
density, but such an operator is not known in general. 
Therefore, we assume that the local formulation of the 
first law of thermodynamics (1231) also holds for the corre- 
sponding operators pi, [H, [26l - l28l |. This leads to the equa- 

tioiQ 



di = Tda 



a)da,i 



(32) 



that expresses the entropy density operator <r(r, t) in 
terms of the internal energy density operator i(r, t) 
and the operators a,i (r , t) for the other relevant vari- 
ables that are known. Notice that the relevant variables 
a; £ {aj, af } are not distinguished into conserved and 
non-conserved variables in this paragraph. Equation (|32[) 
can now be rearranged into 



T 



T 



(33) 



providing a dynamical equation for the entropy density 
operator. Using the definition (|3"Tj) of the time correla- 
tion functions and their symmetry properties described 
in Ref. @, relations between the time-derived correla- 
tion functions CxY(r,r\t,t') with I,F £ {a, e,ai} can 
be derived. A combination of these relations leads to the 
following transport equation for the entropy time correla- 



^(r>-M0 = T(r ^ )T(nt/) 4 g (r>-M0 

c (rr'tf) 
T(r,t)T(r',t>) ai£{ ' > ' > 



(34) 



T{r,t)T{r',t>) 
a\(rM(r',t') 



'—C sai (r,f',t,t') 



T(r,t)T(r',t') (f ' ? ' M) ' 



This equation constitutes the second main result of this 
paper, since it brings about a microscopic expression for 
the entropy time correlation function. 



C. Dissipation 

By the entropy balance equation (|17l) with the quasi- 
current ([2"T]) and the entropy production (|22|) . the dissipa- 
tion function (|29j) is directly associated with dissipation. 
Also the time-derived correlation functions in Eq. (|3"4")l 
can be related to dissipation, when the validity of the 
fluctuation-dissipation theorem [H, [29[ is assumed. For 
systems sufficiently close to thermodynamic equilibrium, 
the fluctuation-dissipation theorem implies the relation 

C a ^ (r, r', t,t') = ^ Xij (f, r\ t, t') (35) 

or equivalently in Fourier space (see appendix I A ip 

Xdj (r, r',u}) = -ui C aiaj (r , r ' u) (36) 

between the time correlation function C aiaj (r,r',t,t') 
and the absorptive response function^ fll| 



4-(f, ?, t, t') = ^([&i(r, t),&j(f, *')])eq (37) 



for any variables di(r, i) and a,j(r,t) including also the 
internal energy density e(r, t) and the entropy density 
a(f,t). In the following, the meaning of the dissipa- 
tion functional *H and the absorptive response function 
Xij{r,f',t,t') in the context of dissipation is considered 
in more detail. 



1. N on- equilibrium dynamics 

The total work W(t) done on a system with the effec- 
tive Hamiltonian H e s(t) is given by [1, H(| 



W(t)=Tr(p(t)H oS (t)) 



(38) 



10 Although Eq. (132 1 provides an expression for the entropy oper- 
ator 5-(r, t), this variable cannot be incorporated into EDDFT 
in the usual way, since the MZFT requires a set of independent 
relevant variables Hill. 



11 Sec Rcfs. 8, 26, 30] for the symmetry properties of the absorptive 
response function Xij(f,r',t,t'). 
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Its rate of change W(t) is the energy dissipated per time 
in this system and identical with the dissipation func- 
tional @, m 



m = W(t) = Tr(p(t)H eS (t)) 

n „ 

= -V A z r ai {r,t)h\{r,t) , 



(39) 



where the averaged relevant variables a,(r, t) now include 
also the internal energy density e(r, t), but not the de- 
pendent entropy density <r(r,t). Notice that Eq. (|3"9")l 
can be understood as a generalization of Eq. (33) in Ref. 
[3. The time integral of Eq. (J2U) is the total amount of 
dissipated energy 

W diss = /dt <K = V fdt f d 3 r di(r, t)a\{r, t) . (40) 

i=1 JR JR 3 



Equilibrium correlations 

Using the fluctuation-dissipation theorem ([33]) , the dis- 
sipated energy (|4U]) can also be expressed in dependence 
of the time-derived correlation functions C aiaj (r, r ', i, t') 
or equivalently in dependence of the absorptive response 
functions Xij{^i t, t'). For this purpose, the symmet- 
ric equilibrium susceptibility matrix (complex response 
function) lllfl 



Xij(r,f',t,t') 



5a.i(f,t) 



Sa)(r',t') 



(41) 



is used in order to express the averaged relevant variables 
(ii(r, t) = a° q (r) + Aaj(r, t) in terms of the corresponding 
thermodynamic conjugates a\ (r, tV^: 

Aa, (r, i) = V fdt' f d V xa (r, f, t, t')a) (r ', t') . (42) 

JR 3 

Inserting this expression into Eq. (|4"U)) leads to 
Wdiss - V fdt fdt' fd 3 rfd 3 r' 

ij = 1 JR JR JR 3 JK 3 (43) 



In Fourier space, the complex response function (|4ip can 
be decomposed as [1, [ll[ 

Xij {?,r',u) = x'ij (r , r", u ) + i x£ (f , r » (44) 



with the (Fourier transformed) reactive response function 
Xij (?, r', uj) and absorptive response function x'L if, u>\ 



Usins 



this decomposition, the symmetry properties 



I26|,l30j Xji(r',r, -u) = Xy(f,f',w) and -x"i(r',r, -uj) 
Xij (r, f ', uj) of the reactive and absorptive response func- 
tions, respectively, as well as Parseval's theorem (see ap- 
pendix IA 2[) , it can be shown that Eq. (|43| is equivalent 
to 

W diss = ±- V fd 3 rfd 3 r'fdw 

x (r, -u)u)Xij (r, u)a) (r", uj) 



(45) 



In accordance with the second law of thermodynamics, 
the dissipated energy Wdiss must be non-negative in a sta- 
ble system (see Sec. Ill C] ). Since this stability condition 
must also hold for any generalized Hclmholtz free-energy 
functional T and hence for an arbitrary set of thermo- 
dynamic conjugates a\(f, uj), the productF^l uj y'L (r, f', uj) 
has to be positive semi-definite for any u [|. As a conse- 
quence of Eq. (|3"6"j). also C aiaj (f,r',uj) has to be positive 
semi-definite. 



IV. EXTENDED DYNAMICAL DENSITY 
FUNCTIONAL THEORY 

So far, a general set of relevant conserved or non- 
conserved variables ai (f , t) with general transport equa- 
tions © and (pTO]) . respectively, has been considered. In 
this section, we now focus on current EDDFT [ll[ with 
its conserved, real, and independent variables a,i{f,t) = 
a£(r, t) with i G {1, . . . , n} including the internal energy 
density and with the corresponding specified transport 
equations. We at first summarize the EDDFT equations 
and the associated transport equations for time correla- 
tion functions. Afterwards, we discuss the derivation of 
EDDFT from a dissipation functional and consider the 
hydrodynamic limit of EDDFT. These considerations are 
valid for general variables ai(r,t) that have to be con- 
served, real, and independent, as usual in the context 
of EDDFT [ll|. For more particular EDDFT equations 
with the concentrations c; (r , t) and the energy density 
i(r,t) as relevant variables, see Ref. (Tlj . 



A. Non-equilibrium dynamics 

In its currently most general form, the EDDFT equa- 
tions are given by (Tl| 



a,i(r,t) + VfJ ai (f,t) = 



(46) 



In the linear regime near equilibrium, Eq. (142 l i follows directly 
from a functional Taylor expansion of the variable a; (r, t) with 
respect to the thermodynamic conjugates a\ (r, t) about the equi- 
librium state, in which the thermodynamic conjugates vanish. 



Notice that the sign of this expression depends on the definition 
of the Fourier transformation. 
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with the currents J Qi (r , t) = (f, t) + Jg ! (r , t) and 

j£(r,i)=Tr(p(t)>)(r,0)) , (47) 
n p 

Jb (f,t) = -Y, dV ^ *)^*"o} (^ ' ■ (48) 

Here, £)WJ (r, r ' i) denotes a diffusion tensor, for which 
explicit expressions can be found in Ref. [Hj |. We 
note that £)(* J ') (r , r ', i) in Eq. (jl5]l contains all non- 
instantaneous contributions. These arc mainly dissipa- 
tive in nature. For certain systems such as, for example, 
nematic liquid crystals, there are also reversible contri- 
butions even in the hydrodynamic regime [f| H| . 

B. Equilibrium correlations 

In the linear regime near equilibrium, the correspond- 
ing transport equations for the equilibrium time correla- 
tion functions C ciiaj (r, r', tj3 are given by [ll[ 

C aiaj (r, f, t) + V P - jW> (r , r', t) = (49) 

with the total currents 3^ (r, r', t) = J^ j \f,r',t) + 
i^\r,r',t) and their contributions 

n 

^ (r, r", t) = -J2 d V fi W (r, r ") C fcj - (r", r ', t) , (50) 
fc=i 3 

n „ 

J<f (r, r ', t) = -J2 d V r W (r, r ") C kj ( f » f, t) - (51) 

For explicit expressions for Qeq (F, r') and riq (r, f"), 
sec Ref. [HI- 

C. Dissipation 

The reversible currents (|4T|) of EDDFT often vanish. 
This is in general the case, when all relevant variables 
di (r , t) have the same time-reversal behavior Q . Then 
only the dissipative currents (|48|) remain and the com- 
plete EDDFT equations (|46p can be derived from a dis- 
sipation functional that is defined by Eqs. (fT8|) and ([29]) . 
Notice that the internal energy density is treated sepa- 
rately from the other relevant variables in Eq. ([29]) , while 
it is included in the relevant variables cii(r, t) in Eqs. (|46|) - 
(|48|) . For instance, the traditional DDFT of Marconi and 
Tarazona [HO] 

p(r,t) + V r -JP(r,t) = 0, (52) 



Here, we made use of the fact that the equilibrium time correla- 
tion functions C ai aj (r, r', t, t') are translationally invariant with 
respect to time so that t' can be omitted. 



JP (r, t) = -(3D p(r, t)V^ (f, t) (53) 

with the one-particle density p(f, t) and the translational 
short-time diffusion coefficient Dq can be derived for- 
mally from the dissipation functional U\ with the dissipa- 
tion function 

**n(?,t) = \pD p{r,t)(V ?P \r,t)f . (54) 

This reformulation of DDFT in terms of a dissipation 
functional constitutes an alternative representation of 
DDFT besides the usual DDFT equations. Such a 
representation by a dissipation functional is very ad- 
vantageous, when some of the relevant variables (in 
this example the one-particle density p(f,t)) have to 
be parametrized by other order-parameter fields and 
dynamical equations for these parameterizing order- 
parameter fields are needed. The derivation of PFC mod- 
els from DDFT, for example, involves such a paramctriza- 
tion of the one-particle density. While this is pretty 
straightforward for isotropic s yste ms with a purely trans- 
lational particle density [U |32| , the derivation of the 
dynamical equations is getting considerably more com- 
plicated for additional orientational degrees of freedom 
both in two [H, [ljl and three [TtJ spatial dimensions 
(for a recent review see [2(|). The usage of a dissipa- 
tion functional instead of a DDFT equation considerably 
simplifies the situation. 



D. Special cases of the EDDFT equations 

In this subsection, we consider two special cases of 
the EDDFT equations (j4"l))) - (j4"8")) in order to demonstrate 
their applicability and to compare them to equations that 
are known from the literature. These special cases are the 
dynamical equations for non-isothermal colloidal suspen- 
sions and the hydrodynamic limit of the EDDFT equa- 
tions. 



1. Non-isothermal colloidal suspensions 

We consider a non-isothermal colloidal suspension, 
whose state is described by a given generalized Hclmholtz 
free-energy functional F[c, T] in terms of two relevant 
variables, the (conserved) concentration field c(r, t) and 
the (non-conserved) temperature field T(f,t). Such a 
functional can, for example, be obtained from static clas- 
sical density functional theory (DFT) [3~4T - l36l | . when a 
local approximation is applied to the DFT functional 
•F(T, [c]) with the spatially homogeneous temperature 
parameter T. Instead of the temperature field T(f,t), 
also the internal energy density e{r, t) or the entropy den- 
sity <j(f, t) lead - in combination with the concentration 
field c(r, t) - to a set of two relevant variables that are 
appropriate to describe the considered system. Here, we 



8 



choose c(r, t) and a(r, t) as the relevant variables and the 
generalized internal energy functional 

£[c, cr] = T[c, T]+ [ d 3 r T(f, t)a(f, t) (55) 

as the corresponding thermodynamic functional. The 
thermodynamic conjugates of the relevant variables are 
thus defined as &{f,t) = 5£/Sc(f,t) and a^(r,t) — 
8£ / 5a{r 1 t) in this section. Since c(r,t), e{r 1 t), and 
cr(r, t) are even under parity inversion and time rever- 
sal, there are no reversible currents fill ]: J^(r,t) = 
J^(f,t) = J^(r, t) = 0. Hence, the dynamics of such 
a non-isothermal colloidal suspension is described by 

c(f,t) + V r J^(r,t) = 0, (56) 
&{r,t) + V r J£{r,t) = Q a (r*, t) . (57) 



simplify considerably. In particular, the diffusion ten- 
sor £)( l J) (r, r', t) in Eq. (|48p and the transport matrices 



These dynamical equations follow directly from Eqs. (|46|) . 
(fTT)) . and (jnj)- The dissipative currents in Eqs. ([5d]) and 
(j57|) result from Eq. (|48|) and are given by 



(r, i) = - / d 3 r' /3^ cc > (f, f', t) c h (r ', f ) 



f d 3 r' /3£>( Cff ) (r , r i) Vf»<r" (r t) , 



(58) 



(59) 



Using these expressions for the dissipative currents 
J£(r, £) and J^{f, t), an explicit equation for the dissipa- 
tive internal energy current (r , i) can be derived from 
Eq. (|26[) . This explicit equation is given by 



(60) 



Jg(f, t)=- I dV (r, r", £)V^c"(r", t) 

d 3 r' /3L>^ (r, r', ^V^a^f', t) 

l 3 

with the diffusion tensors 

D (£c) (r,r',t) = cV,t)£> (cc) (r,r',t) 

+ T(r,i) J D^)(r,f',i), 

J D< £,T ) (r, r', t) = c" (r, t)D (c ^ (f, f", t) 

+ T(r,t)D(™\f,f',t) . 

2. Hydrodynamic limit 



All expressions presented in this paper so far are be- 
yond the scope of hydrodynamics and applicable also for 
non- hydrodynamic wave vectors k and frequencies uj. In 
the hydrodynamic limit (k — > 0, u> — > 0), the transport 
matrices become constant and the transport equations 



(61) 
(62) 



n ( ei'(r,r') and T ( e ^>(f,r') in Eqs. $5$ and (J5TJ, respec- 
tively, are constant in the hydrodynamic limit (see Ref. 
[13 for details). 



V. CONCLUSIONS AND PERSPECTIVES 

For a general set of independent relevant variables, we 
derived microscopic transport equations for the entropy 
density [see Eqs. JT7J), and ([27])] as well as 

for the entropy time correlation function [see Eq. (|34[) ] in 
the framework of linear irreversible thermodynamics. We 
thus complemented current EDDFT |nj with a balance 



equation for the entropy density that has previously not 
yet been considered in the context of EDDFT as a rele- 
vant variable. We furthermore derived the microscopic 
expression (|29[) defining a dissipation functional, from 
which the dissipative dynamics of EDDFT can be de- 
rived. This reformulation of EDDFT in terms of a dis- 
sipation functional complements the EDDFT equations 
and states that also DDFT [l2l - [l6j can be derived from 
a dissipation functional. Our dissipation functional is 
especially useful, if some of the relevant variables have 
to be parametrized and dynamical equations for the pa- 
rameterizing quantities are needed. This is in p articular 
the case in the derivation of PFC models [l7H2(j | from 
EDDFT. 

Our approach can be generalized to include macro- 
scopic variables as well. These are variables, which relax 
on a sufficiently long, but finite time scale in the long 
wavelength limit. Their combination with the hydrody- 
namic approach has been pioneered by Khalatnikov near 
the A transition in superfluid 4 He [37j. In this case the 
degree of order, the modulus of the macroscopic wave 
function, becomes a slow macroscopic variable as this 
second-order phase transition is approached. Later the 
dynamics of the degree of order has also been introduced 
as a macroscopic variable near other phase transitions of 
second or weakly first order. 

Another class of slow variables arises when there is an 
energy scale in the system, which is much smaller than 
all other energy scales. A classical example of this case 
is superfluid 3 He-A, where the magnetic dipolc-dipole in- 
teraction is of order 2 x 10~ 7 K (temperature associated 
with the appropriate energy) and thus very small com- 
pared to all other energy scales including the energy gap. 
As a consequence of the presence of the magnetic dipole- 
dipole interaction the magnetization density is no longer 
a truly conserved hydrodynamic variable, but a rather 
long-lived macroscopic variable, whose dynamics can be 
incorporated into the hydrodynamics of superfluid 3 Hc-A 
[28l |38| , where it turns out to be crucial for the under- 
standing of the NMR spectra in the superfluid phase. 

A class of macroscopic variables, which is particularly 
important for complex fluids including polymers, col- 
loids, liquid crystalline elastomers etc. are strain fields 
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associated with a transient network. In these cases the 
system reacts like a fluid below a certain frequency, the 
so-called Maxwell frequency, and like a solid or glass 
above this frequency. Therefore, one has also explored 
the consequences of a transient strain field on the macro- 
scopic behavior by combining it with the hydrodynamic 
approach [39l - [44j . It therefore appears like a natural next 
step to incorporate macroscopic variables such as a tran- 
sient strain field into the EDDFT approach. 

In contrast to hydrodynamic equations, our results 
are not restricted to the hydrodynamic regime (k — > 0, 
lu — > 0), but they arc also applicable for larger wave 
vectors k and frequencies uj. However, our results that 
base on linear irreversible thermodynamics are restricted 
to systems that are sufficiently close to thermodynamic 
equilibrium, where the local formulation of the first and 
second law of thermodynamics are valid. This is in gen- 
eral the case for all passive systems, but not for active sys- 
tems. The generalization of our results to active systems 
that are far from thermodynamic equilibrium, where the 
entropy production can also be negative locally or in a 
subsystem and linear irreversible thermodynamics is not 
applicable, will therefore be an important task for the fu- 
ture. In this context, especially the crucial question for 
the existence of an entropy density operator will have to 
be answered. 

Entropy has been considered in the context of DDFT 
in recent work of Espahol [5i| and Schmidt J4(| . While 
in our work described in the present paper a generalized 
Hclmholtz free-energy functional is chosen as an appro- 
priate thermodynamic functional, the DDFT equations 
in the work of Espahol base on an entropy functional and 
do not include the entropy density as a thermodynamic 
variable (45[. As a further difference, the balance equa- 
tion for the entropy density and the dissipation functional 
are not used in the approach of Espafiol. Schmidt, on the 
other hand, proposed DDFT equations for a one-particle 
density and an internal energy density on the basis of a 
generalized grand-canonical potential functional that de- 
pends functionally on the one-particle density and the en- 
tropy density and is minimized by both densities in ther- 
modynamic equilibrium [Ifj. The approach of Schmidt 
also includes a dissipation functional from which the dy- 



namical equations for the one-particle density and the 
internal energy density can be derived, but - opposed to 
our approach - its microscopic justification is not empha- 
sized and addressed. 
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Appendix A: Fourier transformation 

Since there are different definitions of the Fourier 
transformation in the literature, here we give the def- 
inition we used in the context of the work presented. 
We further state a useful theorem related to the Fourier 
transformation. 



1. Definition 

Within the definition we used, the Fourier transforma- 
tion of a time-dependent function X(t) is given by 

X(w) = [dtX(t)e iult , 

7 , (AI) 
X(t) = — dco X{u)e~ iut 

with wel. 

2. Parseval's theorem 

If X(t) and Y(t) arc two timc-dcpcndcnt square- 
intcgrable functions, Parseval's theorem states 

fdtX(t)Y(t) = -!- jdujX(uj)Y{uj) , (A2) 

where the bar ~ denotes complex conjugation. 
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